import numpy as np
import pandas as pd
import statsmodels.api as sm

# 创建一些虚拟数据
np.random.seed(1)
x = np.linspace(-5, 5, 100) # 自变量X值域[-5,+5]间均匀分布的数列
y = x + np.random.normal(size=100)

# 添加常量项 (即截距)
X = sm.add_constant(x)

# 使用OLS普通最小二乘法创建并训练模型
model = sm.OLS(y, X).fit()

print(model.summary())